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The dual observables as order parameters for center symmetry are tested at finite isospin chemical 
potential in a Polyakov-loop enhanced chiral model of QCD with physical quark masses. As a 
counterpart of the dressed Polyakov-loop, the first Fonrier moment of pion condensate is introdnced 
for /i/ > m,r/2 under the temporal twisted boundary conditions for quarks. We demonstrate that 
this dual condensate exhibits the similar temperature dependence as the conventional Polyakov-loop. 
We confirm that its rapid increase with T is driven by the evaporating of pion condensation. On 
the other hand, the dressed Polyakov-loop shows abnormal thermal behavior, which even decreases 
with T at low temperatures due to the influence of pion condensate. We thus argue that in QCD 
the critical temperature extracting from a dual observable may have nothing to do with the quark 
confinement-deconfinement transition if the quark mass is very small. 

PACS number(s): 12.38.Aw; 11.30.RD; 12.38.Lg; 


I. INTRODUCTION 

The main phenomena in QCD at finite temperature 
and density are the chiral restoration and deconfining 
phase transitions. In the chiral limit, the standard order 
parameter for chiral transition is the quark condensate. 
However, it is conceptually difficult to define an order pa¬ 
rameter for deconfinement in QCD. Usually, the expec¬ 
tation value of the Polyakov loop (PL) is adopted as the 
indicator for quark deconfining transition. This quantity 
is a true order parameter for center symmetry in pure 
Yang-Mills theory. Nevertheless, this symmetry is badly 
broken by the light dynamical quarks in QCD. These two 
order parameters had been extensively studied in lattice 
QCD. It suggests that both phase transitions are smooth 
crossovers at finite T and zero density and two pseudo 
critical temperatures are very close to each other Si¬ 
lt is well known that the nonzero Dirac zero-mode 
density is responsible for the dynamical chiral symme¬ 
try breaking in QCD, according to the celebrated Banks- 
Casher relation (3| . An interesting question is to what ex¬ 
tent the spectral properties of the Dirac operator contain 
the confinement information. Recently, some authors 
have tried to link the Dirac spectral modes to the PL or 
its equivalent quantities with the same winding number 
in the time direction In these studies, some dual 

observables are introduced as the new order parameters 
for center symmetry by using the twisted boundary con¬ 
ditions for quarks. Especially, it is demonstrated in the 
formalism of lattice QCD [§-[3 that the dressed polyakov- 
loop (DPL) interpolates between the chiral condensate 
and the thin PL. The studies from the functional meth¬ 
ods [HHil and effective models [IM3 also suggest the 
DPL shows the order parameter-like behavior like the 
thin PL. 

In principle, one can construct many dual observables 
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which transform in the same way as the thin PL under the 
center transformation They are true order param¬ 

eters for deconfinement in the static limit m —> 00 , where 
m is the mass of dynamical quarks. However, the cen¬ 
ter symmetry is seriously broken in QCD since the light 
quark masses are very small. Then, a question naturally 
arises: to what extent do these quantities still contain the 
confinement information? Recently, it is demonstrated 
in the NJL model that the rapid rise of DPL near is 
totally driven by the chiral transition [I^. The author 
argues that the reason for it is the lacking confinement of 
NJL. Note that such an explanation is not so convincing 
since the thermal behavior of DPL calculated in NJL is 
quite similar to that obtained in other methods. So even 
the argument is only against the NJL result, the main 
conclusion that the DPL does not reflect the deconfine¬ 
ment transition in [l^ may also be true in QCD. 

To further test weather the dual observables can be 
used as order parameters, we extend the previous study 
to finite isospin chemical potential /i/ by simultane¬ 
ously considering the pion condensation and the twisted 
boundary conditions in this paper. We mainly concen¬ 
trate on the thermal properties of two simple dual ob¬ 
servables for /i/ > TO.n./2 the DPL and the first Fourier 
moment of the generalized pion condensate. Here we re¬ 
fer to the later as the dual pion condensate (DPC), which 
is the counterpart of the DPL at finite /i/. Due to the in¬ 
fluence of pion condensate, the thermal property of DPL 
may change explicitly for > m,r/2. Thus it is interest¬ 
ing to check whether the DPL still behaves like an order 
parameter in this situation. Second, similar to the DPL, 
the DPC transforms in the same manner as the thin PL 
under the center transformation. So it is also interesting 
to explore whether this simple dual quantity can be used 
to indicate the deconfinement transition at finite isospin 


^ It is well known that the charged pion condensation appears at 
low temperatures for fij > m,r/2 in QCD [T^ . 
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density. 

We employ the PL enhanced NJL model (PNJL) in 
our investigation by adoptin g th e U(l)-valued boundary 
conditions. Compared to the advantage of PNJL 

is that the PL dynamics is included to partially mimic 
the confinement. Moreover, the pion condensate and PL 
obtained in lattice simulations |20l.l2ll| for /r/ > to^/2 can 
be well reproduced in this model [22|. The paper is or¬ 
ganized as follows. In Sec.II, the dual pion condensate is 
defined and the PNJL model with the twisted boundary 
conditions for fj,j > to^/ 2 is introduced. The numerical 
results and discussion are given in Sec.III. In Sec.IV, we 
summarize. 


II. DUAL PION CONDENSATE AND PNJL 
MODEL WITH TWISTED BOUNDARY 
CONDITION FOR /x/ > m„/2 


A. Dual pion condensates for /2 

According to the dual quark condensates are de¬ 
fined as 


= ( 1 ) 

where n is an integer and cr((f>) is the generalized quark 
condensate 

a{(j)) = {iPiP)^ =-^{Tr[{m +D^)-^]), (2) 

which is obtained with the twisted boundary condition 
in the time direction 




The in ([5]) is the Dirac operator without the quark 
mass for the twisted angle (j). Note that </> = tt corre¬ 
sponds to the physical boundary condition. The dressed 
PL is defined as the first Fourier moment of cr((^), namely 

(4) 

In the lattice language, this quantity only includes the 
contributions of (infinite) closed loops with the winding 
number one in the compact time direction Q. So it be¬ 
longs to the same class as the thin PL under the center 
transformation. 

In the same way, we can introduce the dual pion con¬ 
densates at finite /x/. The charged pion condensates for 
/x/ > m 7 r /2 are defined as 


=TT+ = -^e 


i9 


{'ipijsT-ip) = tt = —^e' 

v 2 


-ie 


(5) 


where = (ti ± T 2 )IV^ and Ti is the Pauli matrix 
in quark flavor space. In nonzero tt indicates the 


spontaneous breaking of the isospin I 3 symmetry and 
the breaking direction is described by the phase factor 
6 . Without loss of generality, we adopt 0 = 0 in the 
following and the pion condensate is expressed as 

(•0*75X1 x/)) = tt. (6) 

Similar to o, we can define the dual pion condensates 

= (7) 

where 7r(0) is the generalized pion condensate under the 

boundary condition (jS]), which takes the form 

7 ^( 0 ) = -^{Tr[i'y5Ti{m + ( 8 ) 

As mentioned, the DPC is defined as the first Fourier 
moment of 7r{0), namely 

(9) 

Analogous to the DPL (and also the dual density pro¬ 
posed in M), is gauge invariant which merely in¬ 
cludes the contributions of closed loops with wingding 
number one. Thus it is another simple dual observable 
transforming in the same manner as the thin PL under 
the Z(3) center transformation. It is interesting to check 
whether this quantity also exhibits an order parameter¬ 
like behavior with increasing T at finite isospin density. 

The previous studies [2l|, suggest that under the 
physical boundary condition, the pion condensate com¬ 
petes with the quark condensate for /x/ > 77X7r/2. Or in 
other words, the quark condensate partially rotates into 
the pion condensate when the isospin chemical poten¬ 
tial surpasses the half of the pion mass at zero T and 
their competition becomes more involved at finite T. We 
can expect that there may exist the similar interplay be¬ 
tween these two condensates for other twisted boundary 
angles. This implies 7r(0) affects cr(0), and vice versa. So 
the thermal behavior of DPL at jii > m,n./2 may deviate 
significantly from that at zero /x/ due to the influence of 
7 r(0). We will test whether such a deviation still sup¬ 
ports the DPL as an indicator for quark deconfinement 
transition with physical quark masses. 


B. PNJL model for > rriT^ with twisted 
boundary condition 

We adopt the following lagrangian of two-flavor PNJL 
model 


£ = 0(*7^£)7‘-|-7o/x-mo-*A75ri)0 

+9s (00) V (0X75f0)^ -ff^(07/x0)^ 


-9I (007X.0 )^ - U , 


( 10 ) 
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where the last term is the effective PL potential. This 
type model has been used to study the DPL at zero den¬ 
sity M- Compared to , we ignore the eight-quark in¬ 
teraction but include four-quark vector interactions with 
two different couplings. It is demonstrated in [2^ that 
the mismatch between and can lead to non-anomaly 
flavor mixing at finite baryon and isospin densities. 

The rriQ is the matrix of current quark masses 


where Pu{d) is the u(d) quark density. 

In the following, we only consider the situation with 
finite gi and zero /i. In this case, the baryon number 
density is zero for </) = tt and $ equals to $ strictly. Min¬ 
imizing the thermal dynamical potential o, the motion 
equations for the mean fields cr, tt, $ and the density pi 
are determined through the coupled equations 


Too = 


rriy 


rud 


( 11 ) 


dVt 


= 0 , 


dVt 

dir 


= 0 , 


dn 


= 0 , 


dn 

dpi 


= 0 . 


(18) 


and we choose mu = md = m. The /t is the matrix of 
quark chemical potentials 


M = 


Pu 


with 


P = 


Pu + Pd 


Pd 


Pb 


P + Pi 


p- PI 


and 


PI = 


l^d 


( 12 ) 


(13) 


The Pb and pi in (TT^ are the baryon and isospin chemi¬ 
cal potentials, which correspond to the conserved baryon 
and isospin charges, respectively. Following [U, [2^, we 
introduce a small parameter A in (II0|1 . which explicitly 
breaks the symmetry. 

The mean field thermal potential of PNJL model for 
Pi > m^/2 under the physical boundary condition has 
been given in where vector interactions are ignored. 
Its form is slightly modified when considering the vector 
interactions 


fl =W($,$,T) - 2N, 
d^p 


d^p 


[E-+ E;]e{A^ - f) 


-2T 


-l-e 


f d^P r 

J ( 27 r) 3 L 


-3{Ep -E)P 


In 


+ ln 


1 -f 3($ -1- $e"(^»> -^')^)e-(^p 
1 -H 3($ $e“(^p+^')^)e“(^^+''' 


This set of equations is then solved for the fields u, tt, $ 
and Pi as functions of T and pi. 

Under the generalized boundary condition, the mod¬ 
ified quark chemical potential p' in (fHl) should be re¬ 
placed by iT{cj) — tt) which is nothing but an 

effective imaginary chemical potential. Strictly speaking, 
the p' for ^ TT should also contain the density-related 
contribution 2g®(p„ -1- pd) even the real p is zero. This is 
because the imaginary chemical potential also leads to a 
nonzero baryon number density. It has been shown in [lB | 

that the coupling g® only has significant effect on ^ for 
T > l.STc in PNJL. Since we are only interested in the 
thermal behavior of dual observables near and below Tc, 
the correction 2g^{pu + Pd) is ignored in our calculation. 
Note that the p'^ is still real and keeps the form as (1171) . 
The reason is that it is the difference between and 
p'j^ and their imaginary parts cancel each other out for 

According to the definition of DPL(and also the DPC), 
the twisted boundary condition is imposed on the Dirac 
operator and the bracket (• • A still keeps the an- 
tiperiodic condition with A = tt [3, . So in our calcula¬ 

tion, the $ as a function of T and pi is first obtained by 
solving (fT^ using the physical boundary condition. The 
other quantities, such as cf{4>), Tr{<p) p{(j)) are then 
determined by the following coupled equations 


)/5 


dn 

daid) 


= 0 , 


dn 

dTr{(j)) 


= 0 , 


dn 

dpi{(j)) 


= 0 , 


(19) 


+e 3(^;p+A))/3j +in[l-H3($-h$e (^^ (^^ '^consistent with 

+e-3(^^-^')/3l + In [l -f 3($ + $e-(^^+'^')/3)e-(^^+'^')/5 

)'^ I + 5s(o"^ + - 9l(.Pu + Pdf 

~9v iPiJ- ~ Pd) I 

with the quasi particle energy Ep = sj{Ep ± + N'^ 


with $ keeping its value for (j) = tt. Such a treatment is 


C. Model parameters 


and Ep = in which the two energy gaps are 

defined as 


M = m — 2gsU, 
N = X — 2gs'iT. 


(15) 

(16) 


The p' and p'j are the shifted quark and isospin chemical 
potentials 

P' = P- ‘^9viPu + Pd), P'l = PI- ‘^PviPu - Pd), (17) 


(14) calculation, the model parameters related to 

NJL, such as the current quark mass m, the momentum 
cutoff A and the scalar coupling gg are all adopted from 
[^ . which take the values 

TO = 5.5MeV, A = 0.651GeV, and g^ = 5.04GeV■^ 

( 20 ) 

respectively. The vector coupling g"!^ is fixed as 0.255s, 
which is supported by the instanton liquid molecular 
model. 

As for the PL potential, we employ the logarithm form 


24j . It has been reported that this type of U can repro¬ 


duce the LQGD data at finite imaginary chemical poten¬ 
tial, but the polynomial one (which is used in [22|) does 
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FIG. 1: The twisted angle dependences of the quark con¬ 
densate and pion condensate 7r((^) at = 100 MeV for 
different temperatures. 



(b) 

FIG. 2: The twisted angle dependences of the quark con¬ 
densate and pion condensate 7r((^) at = 200 MeV for 
different temperatures. 


not [13 . Following [Tl], the parameter Tq in the loga¬ 
rithm potential is fitted as 200 MeV to reproduce the 
lattice pseudo-critical temperature Tc at zero density. 

For computational convenience, a small A with value 
m/100 is used. It is confirmed that the deviation of our 
main results from zero A is negligible. 


III. NUMERICAL RESULTS AND 
DISCUSSIONS 

A. 0-dependence of quark and pion condensates 

The generalized quark and pion condensates as func¬ 
tions of 0 for /r/ = 100 MeV at different temperatures 
are shown in Fig. [TJ Fig. [TJa indicates that the shapes of 
|cr(0)| for temperatures below and above are quite dif¬ 
ferent: for T = 210 MeV and 240 MeV (or T > T^), the 
quark condensates are concave lines with |(T( 7 r)| < |(t( 0 )|; 
but for T = 150 MeV and 170 MeV (or T < T^), they are 
convex ones with |cr(7r)| > |(t( 0)|. The transition line be¬ 
tween the concave and convex ones takes the wavy shape, 


as displayed in Fig.[T]a for T = 190 MeV (or All 

this is quite different from what obtained in [^, [T 2 , [13 ; 
where only concave curves emerge for vanishing /r and 
/i/. Fig. [Ha also shows that |(t(0 ~ 7r)| first increases 
and then decreases with T due to the impact of pion con¬ 
densate but in SO [13, |(t(0 ^ 7r)| always decreases 
with T. 

In contrast. Fig. |T]b shows that all lines of | 7 r( 0 )| at 
different fixed T are concave curves. We see that in 
the fermionic-like region (namely the area for 0 ^ tt), 
|7 r(0)| decreases with T but it increases with T in the 
bosonic-like region (namely the area for 0 near zero or 
27r). In addition, the curve of |7r(0)| becomes more flat 
with decreasing T. All those is very similar to the 0- 
dependence of the quark condensate obtained for zero /i 
and fii S, O, [13 • The similarity can be understood in 


^ This is also observed in and other chiral model studies. The 
reason for such an anomaly is that the quantity -I- always 
decreases with T but |7r| reduces more quickly since A is zero but 
m is finite. 
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FIG. 3: The temperature dependences of the normalized 
conventional ployakov loop, quark and pion condensates and 
their corresponding dual parters at /ij = 100 MeV. 



T [MeV] 


FIG. 4: The temperature dependences of the normalized 
conventional ployakov loop, quark and pion condensates and 
their corresponding dual parters at = 200 MeV. 


the following way: For /r/ > mT^/2, the quark conden¬ 
sate partially turns into the pion condensate, and thus 
the later inherits some properties of the former. How¬ 
ever, the (/)-dependence of |cr((]ii)| changes obviously due 
to such a transformation, as shown in Fig. [T]a. 

Figure. [5] shows the (/>-dependence of quark and pion 
condensates at different temperatures for nj = 200 MeV. 
Compared with Fig.[Tl we see that is suppressed and 
7 r(0) is enhanced (suppressed) in the fermion-like (boson¬ 
like) region. But the (p and T dependences of these two 
condensates are still qualitatively consistent with that 
displayed in Fig. [TJ 



T [MeV] 


FIG. 5: The temperature dependences of the T-derivatives of 
the conventional ployakov loop, quark and pion condensates 
and their corresponding dual parters at = 100 MeV. 



T [MeV] 


FIG. 6: The temperature dependences of the T-derivatives of 
the conventional ployakov loop, quark and pion condensates 
and their corresponding dual parters at — 200 MeV. 


B. Thermal behaviors of dual condensates 

Two dual condensates and other three (pseudo-) or¬ 
der parameters as functions of T for /r/ = 100 MeV 
and 200 MeV are shown in Figs. ISHU respectively. The 
quark and pion condensates are obtained with physical 
boundary condition cp = ir, which are normalized by 
(To = a(T = 0,/r = 0,/i/ = 0); the dual condensates 
and PL are normalized by their corresponding values at 
T = 240 MeV. We mainly focus on the thermal behaviors 
of these quantities near the phase transitions. 

Figure. 131 shows that \'n{p = 7r)| (<i)) decreases (in- 
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creases) monotonically with T, but \(j{4> = 7 r)| first in¬ 
creases (slowly) up to T ^ 180 MeV and then decreases. 
The similar T-dependences of these quantities are also 
observed in Fig. H) These results are qualitatively agree¬ 
ment with what obtained in by using the polyno¬ 
mial PL potential. As mentioned above, the increase of 
\a{4' = 7 r)| with T is due to the competition between the 
quark and pion condensates. 

Consistent with Fig. [Hb, Fig. [3] indicates that the nor¬ 
malized DPC really behaves like an order parameter for 
center symmetry: analogous to the DPL obtained in 
[l5j |. it keeps rather small value in low tempera¬ 
ture region and gradually becomes larger with T. Like 
the thin PL, the DPC increases monotonically with T. 
However, the normalized DPL in Fig. [3] shows abnor¬ 
mal thermal behavior, which first reduces with T (up 
to T ^ ISOMeV) and then raises. Fig. |3] also shows 
that the DPL even becomes negative near and below 
T ^ 190MeV. The similar T-dependences of DPL and 
DPC are also observed in Fig. 01 

The abnormal thermal behavior of DPL can be traced 
back to the non-concave lines of cr((f>) displayed in 
FigslTDa)-l2Da). As mentioned, due to the influence of 
7 r((/)), |(t((/))| increases with T below (the critical tem¬ 
perature for the I 3 symmetry restoration) for ^ ^ tt. This 
results in the DPL not always raising with T. Actually, 
Figs. [3]0] clearly show that when = 7 r)| increases 
with T, the DPL decreases, and vice versa. So the DPL 
is quite sensitive to the T-dependence of the quark con¬ 
densate. In contrast, Figs. [3]l4] indicate that the DPL is 
insensitive to the thin PL, at least at low temperature 
region. All these suggests that the DPL obtained with 
physical quark masses mainly reflects chiral transition 
rather than deconfinement. Such a conclusion is agree¬ 
ment with the claim given in [l8l| that the rapid change 
of DPL near in NJL is totally driven by the chiral 
restoration. 

Following Ref. [3, we also calculate several suscep¬ 
tibilities which are defined as the T derivatives of the 
quantities displayed in Figs. [3]01 As in [l^, the peak of 
a susceptibility is used to locate the critical temperature. 
The susceptibilities as functions of T for ^7 = 100 MeV 
and 200 MeV are shown in Fig.[S]and Fig. 01 respectively. 

Fig. 0] shows that the PL susceptibility has only one 
peak, which indicates Tf = 163 MeV. But the other 
susceptibilities all have double peaks, one of which coin¬ 
cides with due to the coupling between the PL and 
quark/pion condensate in PNJL. We see that the highest 
peaks of /dT and d\a\/dT in Fig.[5]are very close to 
each other, and the corresponding critical temperatures 
T^"^ and are about 40 MeV lager than . The coin¬ 
cidence of and is consistent with [l^ even the PL 
dynamics and pion condensate are considered in our cal¬ 
culations. In addition. Fig. 0] also shows that the critical 
temperatures T^’^ and (extracting from /dT and 
9 | 7 r|/ 9 T, respectively) almost coincide, which are about 
25 MeV lager than . The slight difference between T'^ 
and can be traced back to a small A = 0.1m used 


in our numerical calculation. All these coincidences are 
also observed in Fig. 0] Actually, the accordance of TJ 
and shown in Figs. 010 ] is completely analogous to 
= T^'^ obtained in [I^l for zero m. We thus get the 
similar conclusion that the rapid change of DPC near 
is driven by the restoration of symmetry, even it ex¬ 
hibits an order parameter-like thermal behavior as the 
PL. 

Note that we also perform the similar calculations by 
varying in this model. We confirm that the thermal 
properties of DPL and DPC shown in Figs. 1-6 do not 
change qualitatively. 


C. The case in the chiral limit 

Here we only show the results with physical quark 
masses. In the chiral limit with finite /i/, the pion con¬ 
densate appears at low temperatures but the quark con¬ 
densate vanishes. Or in other words, the quark conden¬ 
sate totally turns into the pion condensate due to the 
nonzero gi. Correspondingly, the dual quark condensate 
is replaced by the dual pion condensate. In this case, 
the approach used in [T^ for analysing the DPL, such as 
the Ginzburg-Landau method, can be borrowed directly 
to study the DPC. We have checked that without the 
quark condensate, the T-dependence of DPC for finite 
fj,i is much more similar to that of DPL obtained in the 
chiral limit at zero gi [l^. Actually, this case is also 
analogous to the situation with physical quark masses 
and large /i/, in which the pion condensate dominates 
and the suppressed quark condensate can be ignored. 


D. Discussions 

So beyond [l^ , we give further evidences that the dual 
observable may not really reflect deconfinement transi¬ 
tion, even it is constructed from center symmetry. Note 
that in [I^, the author still insists that the DPL calcu¬ 
lated in other methods, such as the formalism of trun¬ 
cated Dyson-Schwinger Equation [T^, can be used as an 
order parameter for deconfinement. However, the dual 
quark condensates obtained in NJL [l^ and DSE [l^ 
are qualitatively consistent with each other. In addition, 
our result and Ref.[l^ all suggest that the T-dependence 
of DPL obtained in PNJL is also quite similar to that cal¬ 
culated in the conventional NJL. So if the rapid change 
of DPL with T merely indicates the chiral transition in 
NJL, there is no reason to explain it as the deconfinement 
transition in other methods. We thus argue that the so 
called coincidence of chiral and deconfinement transitions 
in other formalisms using the DPL as the order parame¬ 
ter may also be an antifact, just as that in NJL El- 

Here we stress that the statement that the DPL only 
reflects the chiral transition given in [l8l| may hold true 
not only in NJL but also in QCD with physical quark 
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masses We argue that the main reason responsible for 
it should be the severe violation of center symmetry due 
to the light dynamical quarks. We point out that it is 
a natural result that the dual observables, such as the 
DPL, contain very limited information on confinement- 
deconfinement transition, unless the dynamical quarks 
are heavy enough. Actually, the abnormal T-dependence 
of DPL displayed in Figs. OH] just reflects that this quan¬ 
tity is a bad order parameter for deconfinement with 
physical quark masses. 

In the heavy quark limit, the DPL approaches the 
PL and both are good order parameters for center 
symmetry. But in the light quark limit, the DPL dif¬ 
fers significantly from the PL. So even the PL has been 
widely adopted to serve as the order parameter for de¬ 
confinement transition, it does not mean that the DPL 
and other dual observables can do this too with physi¬ 
cal quark masses. Actually, even to what extent does the 
PL contain the information of quark deconfinement tran¬ 
sition is also a subtle problem if the quark mass is very 
small. Formally, the definitions of the DPL and DPC in¬ 
volve quark fields which are naturally related to the quark 
and pion condensates, respectively. So it is not strange 
that the rapid change of the former exactly reflects the 
chiral transition in the chiral limit and that of the later 
just indicates the restoration of I 3 symmetry. Thus one 
should be cautious to mention the so called coincidence 
of the chiral and deconfinement transitions through the 
study of dual observables or other quantities constructed 
from center symmetry. 

IV. CONCLUSION 

The dual observables as possible order parameters for 
center symmetry are tested at finite temperature and 
isospin density with physical quark masses. Besides 
the dressed Polyakov-loop, another simple dual conden¬ 
sate, namely the dual pion condensate is proposed for 
/i/ > m7r/2. We investigate the thermal behaviors of 
these two quantities in the PL enhanced NJL model of 
QCD by considering the pion condensation. Our conclu¬ 
sion is that both dual observables contain little informa¬ 
tion on quark deconfinement transition. 

First, we find that the twisted angle dependence of pion 
condensate is quite analogous to that of quark condensate 
obtained at zero ^ and ^/. Correspondingly, the DPC ex¬ 
hibits the similar T-dependence as the conventional PL 
for fij > m,r/2. We demonstrate that the derivative of 
DPC with respect to T peaks exactly at at which 
the pion condensation evaporates. This is very similar 
to the coincidence verified in |1^ . where the critical tem- 

® Actually, a general method, the Ginzburg-Landau analysis is 

used in IlSl to prove the coincidence of the chiral transition and 


perature extracting from the DPL equals exactly to 
in the chiral limit. Thus, we get the analogous conclu¬ 
sion that the rapid change of DPC near is driven by 
the restoration of I 3 symmetry. So even the DPC shows 
order parameter-like behavior, the critical temperature 
extracted from it has nothing to do with the deconfine¬ 
ment transition. 

Second, we find that the DPL displays abnormal ther¬ 
mal property for /i/ > m,r/2, which even decreases with 
T for T < T^. This is quite different from the thin PL, 
which always increases with T. The anomaly arises due 
to the interplay between the quark and pion condensates. 
We verify that the DPL increases with T if the quark 
condensate (its absolute value) decreases, and vice versa. 
This implies the variation of DPL with T is mainly de¬ 
termined by the chiral dynamics rather than the confine¬ 
ment in the situation with physical quark masses, which 
is in agreement with [l8l| . 

We thus conclude that both dual condensates are un¬ 
qualified order parameters for deconfinement transition, 
even they are constructed from center symmetry. Dif¬ 
ferent from [ 1 ^, we argue that this conclusion holds not 
only in the NJL-type model, but also in QCD because 
of the severe violation of center symmetry. We stress 
that if the quark mass is very small, the dual observable 
should naturally loses its role as an effective indicator for 
deconfinement. We thus suspect that the so called coin¬ 
cidence of the chiral and deconfinement phase transitions 
obtained from the thermal properties of dual observables 
may be just an artifact. 

Note that recently the temporally odd lattice has been 
used to investigate the relation between chiral symme¬ 
try breaking and confinement [ 2 ^ [ 2 ^ . It is found that 
the low-lying Dirac zero mode has little contribution to 
the PL, which indicates there is no one-to-one correspon¬ 
dence between the chiral symmetry breaking and confine¬ 
ment in QCD. This conclusion disfavors the prior lattice 
results that the chiral restoration and deconfinement hap¬ 
pen at the same temperature at zero density. So whether 
the chiral and deconfinement transitions coincides or not 
is still under debate and need further study. 

Since there is no sign problem at finite /x/, our study 
can be performed in the lattice calculation. In addition, 
the pion condensation has been investigated in the DSE 
formalism and other effective models of QCD, such 
as the quark-meson model. It is also interesting to in¬ 
vestigate the thermal properties of dual condensates at 
finite /x/ within these methods. 
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